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A key problem in making precise perturbative QCD (pQCD) predictions is how to set the renor¬ 
malization scale of the running coupling unambiguously at each finite order. The elimination of 
the uncertainty in setting the renormalization scale in pQCD will greatly increase the precision of 
collider tests of the Standard Model and the sensitivity to new phenomena. Renormalization group 
invariance requires that predictions for observables must also be independent on the choice of the 
renormalization scheme. The well-known Brodsky-Lepage-Mackenzie (BLM) approach cannot be 
easily extended beyond next-to-next-to-leading order of pQCD. Several suggestions have been pro¬ 
posed to extend the BLM approach to all-orders. In this paper we discuss two distinct methods. 
One is based on the “Principle of Maximum Conformality” (PMC), which provides a systematic 
all-orders method to eliminate the scale- and scheme- ambiguities of pQCD. The PMC extends the 
BLM procedure to all orders using renormalization group methods; as an outcome, it significantly 
improves the pQCD convergence by eliminating renormalon divergences. An alternative method is 
the “sequential extended BLM” (seBLM) approach, which has been primarily designed to improve 
the convergence of pQCD series. The seBLM, as originally proposed, introduces auxiliary fields 
and follows the pattern of the /3o-expansion to fix the renormalization scale. However, the seBLM 
requires a recomputation of pQCD amplitudes including the auxiliary fields; due to the limited 
availability of calculations using these auxiliary fields, the seBLM has only been applied to a few 
processes at low-orders. In order to avoid the complications of adding extra fields, we propose a 
modified version of seBLM which allows us to apply this method to higher-orders. We then perform 
detailed numerical comparisons of the two alternative scale-setting approaches by investigating their 
predictions for the annihilation cross section ratio R e + e - at four-loop order in pQCD. 

PACS numbers: 12.38.Bx, 12.38.Aw, 11.15.Bt 
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I. INTRODUCTION 

Renormalization group (RG) invariance is a central 
principle of quantum field theories: there cannot be any 
renormalization scheme or renormalization scale ambigu¬ 
ity in the predictions for physical observables. However, 
this principle is typically violated if one uses conventional 
scale-setting methods in QCD due to the unavoidable 
truncation of the pQCD series. Thus, a key problem in 
making precise pQCD predictions is how to set the renor¬ 
malization scale of the running coupling properly at each 
perturbative order without introducing any unphysical 
dependence on the choice of renormalization scheme. The 
elimination of the uncertainty in setting the renormaliza¬ 
tion scale in pQCD will greatly increase the precision of 
collider tests of the Standard Model and the sensitivity to 
new phenomena. A review of the QCD renormalization 
scale-setting problem can be found in Ref. ji ij. 

One of the earliest approaches to solve the scale-setting 
problem in pQCD is known today as the BLM ap¬ 
proach, suggested by Brodsky, Lepage and Mackenzie in 
Ref. @]. The BLM approach was inspired by its coun¬ 
terpart in Abelian quantum electrodynamics (QED) 
the Gell Mann-Low scheme Q • In QED, only vacuum- 
polarization insertions contribute to coupling constant 
renormalization. In the BLM work, the optimal scales 
for the running coupling in non-Abelian theory are set 
at each order in analogy to QED by absorbing all of 
the contributions of the /3-function of the process into 
the running coupling. In the BLM procedure, the n/ 
contributions of fermion loops were used to identify the 
leading {/3i}-terms. However this identification is only 
effective up to next-to-leading order (NLO). The BLM 
method was later extended to next-to-next-to-leading or¬ 
der (NNLO) by Brodsky and Lu 3, who also observed 
that the BLM predictions are independent of the choice 
of the scheme if one relates observables to each other 
via commensurate scale relations. However, this method 
cannot be unambiguously extended to pQCD series at 
very high orders, and thus an extra procedure to distin¬ 
guish the {/3,}-terms in a pQCD series is required. 

A proposed extension of BLM is a method called “se¬ 
quential extended BLM” (seBLM) [j|. The purpose of 
seBLM is to optimize the pQCD convergence. In or¬ 
der to accomplish this, a “large /3o-approximation” (§, 73, 
with slight modifications, is adopted to deal with the 
pQCD series. [Here /3o denotes the first perturbative co¬ 
efficient of the running coupling /3-function.] However, 
the seBLM does not distinguish whether the n/-terms 
in a pQCD expansion are in fact related to the renor¬ 
malization group (RG) /3-function, and thus all n/-terms 
are rearranged into {/3i}-terms. Subsequently the {/3j}- 
terms are resummed into the running coupling following 
the scheme of the /3o-expansion. Due to ambiguities in 
determining the coefficients of the {/3i}-terms at higher- 
orders, new colored fields are introduced to fix the se¬ 
BLM scales. A detailed review of the seBLM procedures 
will be given in this paper; we will also point out several 


shortcomings of the seBLM method, which constrains its 
current applicability. 

If one can unambiguously identify all of the {/3j}-terms 
in a pQCD expansion, then they can be systematically 
eliminated by shifting and thus setting the renormal¬ 
ization scale at each order in pQCD. The remaining 
pQCD series will then match the terms of the corre¬ 
sponding “conformal” series with /3 = 0. This is the 
main idea behind the “Principle of Maximum Confor¬ 
mality (PMC)” Plj. The PMC explains why the BLM 
method works so successfully at the NLO level, and it 
can be applied for any high-energy processes up to any 
order [l3|. For some recent N 3 LO and N 4 LO examples 
see Refs, fcil [Til | . 

When one applies PMC scale-setting, the scales of the 
running coupling in the pQCD series are shifted such that 
all contributions related to the /3-function are resummed 
and only the coefficients which are RG invariant remain. 
These coefficients are called the “conformal series” since 
they are identical to the coefficients obtained if the /? 
function were zero. The resulting PMC predictions are 
thus scheme-independent - predictions under different 
schemes only differ by the appropriate shift of the renor¬ 
malization scale. Thus the PMC numerical prediction is 
insensitive to the scheme choice. The PMC predictions 
also have the property that any residual scale dependence 
on the choice of (initial) renormalization scale is highly 
suppressed, even for low-order predictions. The scheme- 
independence of the PMC predictions is also confirmed 
by commensurate scale relations. An all-order demon¬ 
stration of commensurate scale relations in pQCD can 
be found in Refs. 00- The PMC thus obeys stan¬ 
dard RG-invariance and satisfies all RG-properties |20| . 
Furthermore, the pQCD convergence of pQCD series is 
improved due to the elimination of divergent 
renormalon terms. 

Two ways to identify the nonconformal {/3j}-terms 
and to implement the PMC have been suggested. One 
method is based on the PMC-BLM correspondence 
(PMC-I) @,0 , and the other more recent method is a 
theoretical improvement based on the R^-scheme (PMC- 
II) 0 , 0 |. Both of these methods satisfy RG-invariance, 
but they resum the pQCD series in different steps. We 
have shown that PMC-I and PMC-II are numerically 
equivalent by comparing several high-energy processes up 
to the four-loop level 0, [HI . A systematic demonstra¬ 
tion of this equivalence will appear soon [22|. In the 
following, we will adopt PMC-II for the discussions and 
simply call it PMC. 

Recently, an approximate comparison of PMC and se¬ 
BLM was published by Kataev and Mikhailov [23); they 
concluded that the PMC method is “questionable”. How¬ 
ever, we have found that their conclusions were based on 
a misunderstanding and a misuse of the PMC. It is clearly 
important to clarify these issues so that any misuse of 
the PMC method will be avoided. In this paper, we will 
analyze R e + e - up to four-loop level, making a detailed 
comparison of the PMC and seBLM predictions. Since 
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the seBLM is currently not applicable at the N 3 LO level 
and beyond, we will suggest a new method for extending 
it to higher orders, based on lessons from the PMC. We 
will refer to this modified seBLM as ‘MseBLM’. 

The remaining parts of this paper are organized as fol¬ 
lows: In Sec. II, we will present a mini-review of the 
BLM-like scale-setting methods, BLM, PMC and seBLM. 
The main features of these BLM-like methods will be dis¬ 
cussed, and an overview of the differences between the 
PMC and seBLM procedures will then be presented. We 
will show that the PMC and seBLM methods have some 
common features; in particular, the PMC and the seBLM 
/3-patterns are exactly the same. We will also present 
MseBLM in this section. In Sec. Ill, we will present 
numerical results for R e + e - up to the four-loop level to 
provide a detailed comparison of the PMC, seBLM and 
MseBLM predictions. Sec. IV is reserved for a summary 
and conclusions. Two appendices provide further com¬ 
putational details for seBLM. 

II. BLM-LIKE SCALE SETTINGS 

Consider the nth-order pQCD approximant g n for a 
typical physical observable g, 

n 

Qn = r 0 a p (g) + Y n(g)a p+l (g), ( 1 ) 

i=1 

where a = a s /4-7r and g stands for the initial choice of 
scale. The coefficient ro is the tree-level result, p is the 
power of the coupling associated with the tree-level term, 
and r, (g) is the coefficient of the Ah-loop correction. 

In conventional scale-setting treatments, the renormal¬ 
ization scale is fixed to an initial guessed value which is 
usually taken as a typical momentum flow of the process. 
If one uses this method, the prediction is scheme depen¬ 
dent. Furthermore, the scheme and scale dependence of 
the coefficient r^g) and the coupling constant do not 
exactly cancel at any fixed-order, leading to well-known 
ambiguities. In contrast, the BLM, PMC and seBLM 
scale-setting methods improve the pQCD predictions by 
eliminating such scheme and scale ambiguities. 

A. Basic arguments of BLM 

The renormalization scale for the running coupling is 
unambiguously set in QED by summing all vacuum po¬ 
larization contributions, both proper and improper, into 
the photon propagator. Thus a(t) = in the Gell 

Mann-Low scheme Q, where t is the virtuality of the 
exchanged photon and 11(f) sums all n/ fermion loop 
contributions. The running of the QED effective cou¬ 
pling is due to vacuum polarization alone - only vacuum- 
polarization insertions contribute to the effective cou¬ 
pling i- Following this observation, Brodsky, Lepage 
and Mackenzie pointed out that the n/-dependence of 


the pQCD series can also be used at low orders as a guide 
to identify the /3 0 and j3\ terms and thus set the scale of 
the pQCD prediction up to N 2 LO [2|. This simple and 
straightforward n/ method can be applied to processes 
that do not involve the three- or four-gluon couplings at 
leading order. The BLM method also ensures that the 
pQCD predictions analytically continue at N c —> 0 cor¬ 
rectly to Abelian theory [24[ ■ 

We will take the NLO pQCD prediction for a typical 
single-scale physical observable g\ to illustrate the BLM 
procedure. At NLO level, the pQCD prediction can be 
re-expanded as 

gi = r 0 a p (g ) [1 + ( An f + B) a(g)\ . (2) 

The ny-term is due to light quark vacuum-polarization 
insertions. Any initial scale choice and any renormaliza¬ 
tion scheme, including dimensional regularization, can 
be used for the prediction. The coefficients A and B are 
in general different under different schemes; however, the 
BLM prediction g\ is unchanged for any choice of scheme 
due to commensurate scale relations [41] . 

All of the n /-terms can be resummed into the running 
coupling with the help of the one-loop a s -running cou¬ 
pling 

a(/‘BLM) i + /3 0 a{g) In (gl LM /g 2 )' ^ 

where /So = 11 — 2n//3. One then obtains 

pi = r 0 a p (p B LM) [1 + ria(pBLM)], (4) 

where pblm = gexp(3A/p) and r\ = ^A + B. The 
BLM scale /tblm is thus determined solely by A. The 
term 33A/2 in r* serves to remove those contributions 
which renormalize the running coupling, and the result¬ 
ing o (/ xblm ) is the predicted value of the running cou¬ 
pling. Eq.(|3]) indicates that the BLM coupling ci (/ xblm ), 
and thus the BLM prediction g i, is independent of the 
initial scale g , as is readily checked. 

This approach of using the n/-terms as a guide to re¬ 
sum the series through the RG-equation of a s cannot be 
unambiguously extended to higher orders. One reason 
is that the n/-series and the {/3j}-series are not a pri¬ 
ori one-to-one. Another issue is that n /-terms appear at 
higher orders from loops which are ultraviolet finite but 
are not associated with the /3-function of the running 
coupling. Thus reactions with multi-gluon couplings are 
more difficult to analyze using BLM because quark loops 
appear in high-order corrections to the multi-gluon ver¬ 
tex as well as in the propagator insertions [25(. Scale 
setting for the BFKL Pomeron intercept provides such 
an example [13, [ItJ ■ 

Thus, it is necessary to modify BLM at higher-orders. 
We will discuss two suggestions, PMC and seBLM. We 
shall first present an overview of those two suggestions, 
and then present the PMC and seBLM features and their 
consequences sequentially. In this discussion we will ig¬ 
nore quark mass terms and their renormalization. 
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B. An overview of PMC and seBLM 

The purposes of PMC and seBLM are different. The 
PMC is designed to solve the renormalization scheme- 
and scale- ambiguities, whereas the seBLM is designed to 
improve the pQCD convergence. Both the PMC and se¬ 
BLM utilize the {/3j}-series to achieve these goals, rather 
than the simpler n/-terms. In the case of processes where 
the {/3j}-terms for the quark anomalous dimension and 
the QCD /3-function are entangled with each other, ex¬ 
tra steps have to be taken to distinguish those {/3 i}- 
terms fl5l ]. 

When one applies the PMC or seBLM, two steps are 
needed to fix the renormalization scale. The first step in 
both cases is to fix the /3-pattern at each perturbative or¬ 
der and determine the coefficients of all {/3j}-terms in the 
/3-pattern. However, PMC and seBLM use quite different 
methods to accomplish this step: 

• The PMC observes that the /3-pattern at each or¬ 
der originates from a specific pattern and superpo¬ 
sition of the {/3j}-terms coming from all the lower- 
order a s -factors, due to its running behavior. The 
running behaviors are governed by the fundamen¬ 
tal RG-equation, and one can thus identify the /3- 
pattern up to all orders without any ambiguities. 
As we shall discuss below, this procedure can be 
carried out systematically by generalizing the defi¬ 
nition of MS dimensional regularization to include 
extra subtraction terms <5^=1 2 ,—■ The coefficients 
of 2 ,'-. then isolate the {/3j}-terms to a partic¬ 


ular a s -order 00. 

The coefficients of the {/3;}-terms at each order 
can also be determined from the n/-power series 
at the same order, which are calculated under a 
certain renormalization scheme, such as the MS- 
scheme The expressions for the /3o,i,2,3 as 

a function of nj in the MS-scheme can be found 
in Refs. j29l-l37j. There is a subtlety regarding n/- 
terms that are unrelated to the a s -renormalization, 
which must be kept unchanged when applying the 
PMC. Special degeneracy relations among different 
order terms ensure the exact one-to-one correspon¬ 
dence between the n/-terms and the {/3j}-terms at 
the same order, so that all {/3j}-coef5cients can 
be unambiguously fixed up to all orders [HI- It 
has been demonstrated that the degeneracy rela¬ 
tions are not specific to dimensional regularization 
schemes, but are general features of perturbation 
theory M |. 

• The seBLM fixes the /3-pattern by identifying the 
equivalent /3o-powers of the {/3j}-terms via the re¬ 
lation /3 i ~ /3q + 1 - All of the possible {/3j}-terms, 
whose equivalent /3o-powers are equal or less than 
the maximum /3o-power of the considered order, 
form the /3-pattern at each order. In fact the se¬ 
BLM /3-pattern is the same as the PMC one. Ac¬ 
cording to the RG-equation, the running coupling 
at different scales are related by the following dis¬ 
placement equation: 


i(/x 2 ) = a(/xi) - /3 0 ln ( ) a 2 (/zi) 


/3q In 2 


—A) ln d 


1/30/3, In 2 (j|)-/3 2 ln 




1 


\A) j 


J (Ml) + 


\Aj j 


a 4 (/zi) + 0 (a 5 ), 


(5) 


r 


where /z, and /z 2 are two arbitrary scales. At each 
order, the equivalent /3o-powers are the same for 
all {/3j}-terms. For example, at order a 4 , /3 2 ~ 
/3o/3i ~ /3q . Because the maximum /3o-power at 
each order is the same for PMC and seBLM, the 
superposition of the a^-displacement from all lower 
orders, which is used in PMC, will, as we will see, 
result in the /3-pattern of seBLM. From this point 
of view, the PMC provides the underlying principle 
for the seBLM /3-pattern. 

The seBLM coefficients of the {/3i}-terms are also 
determined from the known ny-power series, but in 
a quite different way than the PMC: At the NLO 
level, only the /3o-term needs to be determined, and 
it can be directly fixed by the rij-term. At the 


NNLO level, the /3g-term can be fixed by the nj- 
term, but the /3,- and /3o-terms cannot be unam¬ 
biguously fixed by the n^-term alone, since both 
/3o and /3, are linear functions of n j. To solve this 
problem, the seBLM method introduces n g mul- 
tiplets of fermions in the adjoint representation of 
the color-group, resembling the gluino of supersym¬ 
metric Yang-Mills theory (which makes part of the 
minimal supersymmetric standard model, MSSM). 
The coefficients for the /3,- and /3o- terms are thus 
fixed by recalculating pQCD series with h g gluinos 
and by using the resulting n j- and n g - terms. The 
gluinos are introduced as a technical device to fix 
the {/3i}-coefficients; however this procedure intro¬ 
duces process-dependent complex calculations into 
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higher-order QCD corrections. Thus in order to ap¬ 
ply the seBLM, the pQCD corrections need to be 
recomputed with the new fields in order to extract 
the h ff -dependence. At present, the n ff -dependent 
/3j-function is known to three-loops; i.e ., up to fa, 
and only the Adler D-function has been calculated 
with n g -dependence up to NNLO level. Thus se¬ 
BLM can only be applied at this time up to NNLO 
for D-function-derived processes, such as R e + e - 
and the Bjorken polarized sum rule [23j . 

For the second step, one needs to find the correct way 
to resum the relevant {/3j}-terms, determined from the 
first step, into the running coupling. The PMC and the 
seBLM take quite different paths: 

• The PMC notes that only those {/3j}-terms which 
are related with the renormalization of the run¬ 
ning coupling should be absorbed into the running 
coupling. They are eliminated through the RG- 
equation, and the resultant PMC scales are func¬ 
tions of the running coupling and are in general 
different for different orders. One can choose any 
initial renormalization scale to do the pQCD calcu¬ 
lation as long as its value is large enough to ensure 
pQCD applicability. The final resummed result, 
however, has no or very small residual dependence 
on the choice of the initial renormalization scales. 

• The seBLM treats all the {/3j}-terms on equal foot¬ 
ing without distinguishing whether they should be 
incorporated into the renormalization of the run¬ 
ning coupling. The seBLM scales are obtained by 
adopting the “large /3o-approximation” [i| Q; i.e. 
using the equivalent /3o-powers for the {/3i}-terms 
as a guide. The highest /3o-powers are eliminated 
first, then the one-order-lower /3o-terms, and so on. 
Thus the seBLM scales are effective scales, and one 
cannot a priori decide whether they will lead to 
the optimal behavior of the running coupling. In 
fact, one finds that the seBLM convergence is not 


as expected, even for the quantity R e + e ~. Thus an 
extra treatment, called x-BLM, to further improve 
the seBLM pQCD convergence was suggested in the 
original seBLM paper [j|. 


C. PMC scale setting 

PMC scale-setting provides the underlying principle 
for BLM, a procedure that can be unambiguously ap¬ 
plied to any order HIM El- The PMC utilizes the iden¬ 
tified {/3j}-terms and the RG-equation to determine the 
value of the running coupling at each order and thus the 
“physical” scales of the process. At the NLO level, one 
only requires the /3o-term; absorbing the /3o-term into the 
scale of the running coupling is equivalent to eliminate 
the n^-term, which explains why BLM works so well at 
the NLO level. 

A key step is to identify the /3-terms in the pQCD 
prediction, thus distinguishing the “nonconformal” ver¬ 
sus the “conformal” ^-independent terms. To do this, 
the PMC introduces a generalized dimensional renor¬ 
malization scheme, the /^-scheme, where an arbitrary 
constant —5 is subtracted in addition to the standard 
subtraction (ln 47 r — 7 e) used in the MS-scheme. The 
^-subtraction defines an infinite set of MS-like schemes. 
The PMC scales for different Ra-schemes, e.g. Rg x - and 
i?a 2 -schemes, differ only by a factor e^ 5l_52 ^ 2 . This scale 
relation ensures the scheme independence of pQCD pre¬ 
dictions among different schemes. Moreover, the scale 
displacement between couplings in any i?a-scheme reveals 
all the {/3i}-terms pertaining to a specific order mill. 
By collecting all of the {/?j}-terms that occur at a given 
order, one can identify the /3-pattern of the RG, and thus 
compute each PMC scale order by order, 

According to PMC scale-setting, the explicit /3-pattern 
at each order for the pQCD prediction of the observable 
g can be written as 


Bn = r 0 ,o + rifia(p) + [r 2 , 0 + far 2 ,i] a 2 (/r) + [ 7 - 3,0 + far 2 ,i + 2/3 0 r 3 ,i + /3o7" 3 , 2 ] a 3 (/z) 


+ 


r 4 ,o + far 2 ,i + 2/3i 7 - 3,1 + -/3i/3 0 7- 3 ,2 + 3/3 0 r 4 ,i + 3/3gr 4 , 2 + /3or 4 , 3 




( 6 ) 


where the 770 are the conformal coefficients, while 0 
are the ones related to the running coupling renormaliza¬ 
tion. The degeneracy relations among the non-conformal 
coefficients at different orders are implicitly adopted. 
The /3-pattern fixed by Eq.j 6 | at each order is dictated 
by the RG-equation. It is thus natural to call such a 
/3-pattern the RG-/3-pattern. It is noted that the RG-/3- 
pattern is determined from a superposition of displace¬ 


ments in the running couplings at each order; thus only 
those {/3i}-terms which determine the correct running be¬ 
havior of the coupling are kept. There are cases where 
the coefficients r,;.j are exactly zero for the {/3j}-terms; 
i.e. there are “missing” {/3i}-terms in specific processes. 
This only indicates that those terms have no contribu¬ 
tions to the renormalization of the running coupling. 
This explains why in Ref.@], what the authors refer to 
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as the “correct PMC”, is actually an invalid procedure. 
Since the authors [23 ] transform all rif -terms into {A}~ 
terms which brings unrelated {/3j}-terms into the RG-/3- 
pattern, thus explicitly breaking the PMC procedure. 


Using PMC scale setting, which follows the pattern 
dictated by the RG, all lower-order running couplings 
can be resummed into effective couplings a k (Q k ) as 


where /3(a) 


i 

/U 


Rk, 


Aj^a) 


r k ,oa k (g) + r k ,o k a k 1 (fj,)P(a) ^R k ,i + A^\a)R k ,2 -t-lA[ n 1 \a)R k , n j = r kfi a k (Q k ) , 

= — a 2 (/r) o A a *(M) an( l Qfc is the PMC scale at a fc -order, given by 

Rk ,i + Aj/ "* (a)R kt 2 + ■* ( a)R ki z + • • • 

1 + + (Ai 1} (a)) 2 (i? fe , 2 - R^) + A^a)^ + • • • ’ 

= (—1 )i r -h±M. 


1 

2 ! 


r k , o 


AfV) = ! 


9 2 /3 (d/3 


da 2 


9a 


P <9/3 


,P 2 


P~aZo. + AT + 3(fc - 1)-^- + (fc ~ 1)(* - 2) — 


a 9a 


(7) 


( 8 ) 

(9) 

( 10 ) 


This shows that the PMC scales are in general different 
at different orders, as in QED. This is clear since the /?- 
terms which control the behavior of the running coupling 
and the physical scales which set the virtuality of the 
propagators, as well as the number of effective flavors 
n/, are usually different at each order. 

The final pQCD PMC prediction for g then reads 

Qn = r 0 , 0 +ri, 0 o(Qi)+ r 2 , 0 a 2 (Q2) 

+^3,0O 3 (£?3) + 1"4,oa i (Q4 : ) + • • • . (11) 

At the four-loop level, Q 4 remains unknown, since we 
need to know the five-loop coefficient r' 5.1 to fix its value. 
This leads to some minor residual scale dependence, sup¬ 
pressed by the highest power in a. A practical choice of 
Q 4 is Q 3 . As shown by Eq.©, the PMC scales them¬ 
selves are perturbative series, which also introduce resid¬ 
ual scale-dependence; however, this dependence is highly 


suppressed, even for lower-order analyses. Since the non- 
conformal contributions are absorbed into the scales, the 
PMC predictions have optimal pQCD convergence due to 
the elimination of divergent renormalon terms. More de¬ 
tails and the properties of PMC may be found in Refs. Il¬ 
ia and in a recent short review in Ref. [2ll] . 


D. seBLM and its modified version MseBLM 

As mentioned in the introduction, seBLM can only be 
applied at present to Adler D-function-derived processes 
up to three-loops 0, HTsI frit . In this subsection, we will 
take the three-loop pQCD prediction of a physical ob¬ 
servable g to illustrate the seBLM procedure. Following 
Eq.Jl]), the pQCD prediction g n for g up to three-loop 
level can be written as 


Qn = r 0 + n (a(n) + r 2 a 2 (fj.) + r 3 a 3 (g.) - ) 

= r 0 + n (a(/z) + (A) • r 2 [ 1 ] + r 2 [ 0 ]) a 2 (/z) + (/3 q 


where ro is free of strong interactions, and for conve¬ 
nience, an overall factor r 1 is factored out of the pQCD 
corrections. The /3-pattern at each order in the sec¬ 
ond line is suggested in Ref. @. The first argument no 
of the coefficients r n [n 0 , tii, • • ■ ] correspond to the Ar 
power, whereas the second one ni corresponds to the 
/3i-power, etc. We have omitted the arguments of the 


r 3 [2] + Al • r 3 [0,1] +/3 0 • r 3 [l] +r 3 [0]) a 3 (/z) H-) , (12) 


{A}-coefhcients for brevity. In order to shorten the no¬ 
tation even further: if all the arguments of the coeffi¬ 
cient r n [- ■ ■ , to, 0 , • • • , 0 ] to the right of the index to are 
equal to zero, then we will omit those zero arguments for 
simplicity and write instead simply r n [■ ■ ■ , to] . The coef¬ 
ficient r n [no, ni, ■ ■ ■ ] is usually scale-dependent through 
terms dependent on In Q 2 /fj, 2 , where Q stands for the 
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typical momentum flow of the process. 

The seBLM /3-pattern, shown in Eg. (fl2l) . is fixed by 
using the relation, /?* ~ /3q +1 . By counting the equivalent 
/3o-powers for all the possible {/3j}-terms, the /3-pattern is 
determined by following the decrement of the equivalent 
/3o-powers. For example: At the N 2 LO level, we have 
only /3q ; at the N 3 LO level, we have /3q, /3i ~ /3q and /3o; 
at the N 4 LO level, we have /3q, P2 ~ Po> P1P0 ~ Pq, Pq, 
Pi ~ /3g, /3o; at the N 5 LO level, we have Pq, /3 3 ~ Pq, 

/3 2 /3o - /3 0 4 , A 2 - Pot P 1 P 0 ~ Ao> /3 o, P 2 ~ /3 0 3 , ft A) ~ /3 0 3 , 
/?o> Pi ~ /?o> An etc. 



FIG. 1: Diagrammatic illustration for the ^-pattern of the 
Os-correction contributing to the photon polarization oper¬ 
ator II [§)], following the relation /3i ~ /3J +1 : (a) diagrams 
without any intrinsic renormalization contributions to r 3 [0]; 
(b) diagrams contributing to the /3o-term 7 * 3 [1]; (c) diagrams 
contributing to the do-term 73 ( 2 ]; (d) diagrams contributing 
to the /3i-term 7-3 [0, 1]. The ellipses denote other diagrams 
with the gluon/vertex renormalization. 

A diagrammatic illustration of the seBLM /3-pattern 
can be motivated by the a s -corrections to the photon 
polarization operator n p|. It is based on the fact that 
the ^-functions emerge together with certain a s -powers; 
i.e., (Pqcl), ( P\a 2 ), , (/3,a l+1 ); at a fixed a s -order, all 
possible combinations of /3j-functions can thus be enu¬ 
merated. Fig. m shows the a 3 -corrections to the gluon 
polarization operator II, in which the diagrams for the co¬ 
efficients 9-3 [0], r 3 [1], r 3 [2] and r 3 [0,1] are presented. An¬ 
other explanation of the seBLM /3-pattern can be based 
on an analysis of the Casimir structures of the Adler D- 
function; see Ref. (23|. Those two explanations for the 
seBLM /3-pattern are illuminating; however, they cannot 
explain why such /3-pattern would be applicable to all 


high-energy processes. Note that the seBLM /3-pattern 
is the same as the PMC RG-/3-pattern. As has been 
discussed in Sec lIIBl the PMC provides the underlying 
reason for the seBLM /3-pattern. 


1. Determination of seBLM {ft^-coefficients 

In practice, the seBLM {/3i}-coefficients r n {n 0 , m, ■ ■ ■ ] 
are fixed by using the full n /-power series at the same or¬ 
der without distinguishing whether those n/-terms per¬ 
tain to the RG /3-function. As we have noted, this is the 
underlying reason why the seBLM differs from the PMC. 

At the two-loop level, the /3o-coefficient r 2 [l] can be 
fixed by the linear n/-term C 2 ,in/ in the two-loop n/- 
power series 

r 2 = C 2 ,in/ + c 2 ,0 , (13) 

where the expansion coefficients C 2 ,o and c 2 1 for the case 
of Adler H-function can be found in Ref. [41] . 

At the three-loop level, the /3 2 -coefficient r 3 [2] can be 
fixed by the squared n 2 -term c 3j 2 ri/; however, the Pq and 
Pi coefficients r 3 [l] and r 3 [0,1] cannot be unambiguously 
fixed by a single linear n/-term c 3j in/ in the third-order 
n/-power series 

D 3 = c 3j2 n/ + c 3j in/ + c 3j0 , (14) 

where the expansion coefficients c 3j o, c 3j i and c 3 ? for the 
case of Adler H-function can be found in Ref. (41j. 

To solve this problem, the seBLM method requires the 
recalculation of the process with the help of multiplets 
of strongly interacting MSSM gluinos. The purpose of 
introducing such MSSM gluinos is to determine the {/3j}- 
coefficients from the n/- and h g - power series; in the 
final step, one can set h g to zero in order to obtain the 
numerical predictions of seBLM for pQCD. 

By recalculating the pQCD series with n g gluinos, the 
third-order coefficient changes to 


h = C 3 , 5 n g nf + c 3j4 n 2 + c 3j3 % + c 3j2 n/ + c 3 pn/ + c 3i0 . 


(15) 


The coefficients po and Pi with MSSM gluinos are linear 
functions of nj- and hg [ 43 , thus the {/3o,/3i}-coefficients 
r 3 [l] and r 3 [0,1] can be fixed by c 3)3 and c 3) i. Note 
that the /3-series has only four independent coefficients, 
whereas the n / and h g series given above have six expan¬ 
sion coefficients. Thus not all of the Cy coefficients are 
independent. We can in fact find two relations among 
c 3j 5 , c 3> 4 and c 3i2 ; i.e. 

c 3,5 = 4 c 3,2 • c 3 ,4 and 4 T^c 3i4 = C^c 3i2 , 


which agree with the analytic results derived in Ref. j4l| . 
Thus, the /3Q-coefficients r 3 [2] can be fixed by c 3j s or c 3j4 
or by c 3j2 . 


These seBLM procedures can be extended to all orders. 
For example, the seBLM /3-pattern at the four-loop level 
can be written as 



















Po ' ^4[3] + PoP\ ■ 1] + /3 2 • r 4 [0,0,1 ] +/3q • r 4 [2] + Pi • r 4 [0,1] + /3 0 • »* 4 [1] + r 4 [0]. 


(16) 


The pQCD n/-power series at the four-loop level can be 
written as 

r 4 = c 4j3 n 3 + c^ 2 n * 2 + c 4 pn/ + c 4j0 . (17) 


Again it cannot fix all the {/3j}-coefficients. By recal¬ 
culating r 4 with fig gluinos, the fourth-order coefficient 
changes to 


J 


r 4 = c^ghgrij + C 4 )S ,h 2 g rif + c^figUf + c^^n 3 + c 4i5 fig + c 4i4 n g + c 4i3 n 3 + c 4i2 n 3 + c 4 , 4 n/ + c 4 , 0 - (18) 


The {/3o,/3 4 (--coefficients r 4 [l] and r 4 [0,1] can be fixed 
by c 4j4 and c 4 p; the {/3 0 /3i,/3 2 , ^(--coefficients r 4 [l,l], 
r 4 [0, 0,1] and r 4 [2] can be fixed by c 4j 7 , c 4j 5 and c 4j2 . The 
/9-series msi) has seven independent coefficients, whereas 
the rif and n g series (1181) have ten expansion coefficients. 
Thus there are three constraints which constrain the Cij 
coefficients at this level. Using the equivalence of (flbl) 
and m, we obtain 

9c 4j 3 • C 4j 6 = c 4j g • C 4j g, 

Ca c 4,3 = 87’))c 4 .6 , 

27kc 4) 8 = CUc 4) 9 . 

Thus, only one of c 4j g, c 4i g, c 4i 6 and c 4j3 is independent; 
it can be used to fix the /3g-coefficient r 4 [3]. At present, 
the fourth-order coefficient f 4 is not available for any ob¬ 
servable; the general relations given above can be used 
as a cross-check on 7 4 when it is calculated in the future 
and as a self-consistency check on the seBLM idea. 

The hg-related coefficients are process-dependent. One 
requires the h s -dependent calculation in order to fix all 
the {/3j}-coefficients. This unavoidably introduces extra 


loop calculations into the already complex higher-order 
QCD corrections. Due to the present limited knowl¬ 
edge of the fig-dependent {/3j}-expression and the fig- 
dependent pQCD series, the seBLM can only be applied 
to deal with Adler D-function-derived processes, which 
is now only known at the three-loop level. For future us¬ 
age, we present the seBLM coefficients r'k[l. m, ■ ■ ■ ] at the 
scale p = Q for the Adler D-function in the Appendix [Al 
where Q stands for the measured physical scale. 


2. Determination of the seBLM scales 

After fixing the /3-pattern and the {/3i}-coefRcients for 

a process, the “large /3o-approximation” can be adopted 
with a slight alteration to analyze the pQCD series. This 
first step can be done at any order. As an explanation 
and for the subsequent use of MseBLM, we will present 
a four-loop analysis of the pQCD series. Up to four-loop 
level, the pQCD approximant g n can be formally written 
as 


Qn = r 0 + n 


a{0) + ( Po ■ ?" 2 [1] + r 2 [0]) a 2 (n) + f /3p • r 3 [2] + /3i • r 3 [0,1] + f3 0 ■ r 3 [l] + r 3 [0]) a 3 (/x) 


( pi ■ r 4 [3] + /3 0 /3i • r 4 [l, 1] +/3 2 • r 4 [0,0,1] + /3g • r 4 [2] + /3 4 • r 4 [0,1] + /3 0 • r- 4 [l] +r 4 [0]^) a 4 (/r) H- (19) 


We note that in Ref.[HI], the elimination of {/3i}-terms 
starts with the so-called RG-invariant g n \n=Q- This is 
not a strict condition since the fixed-order pQCD approx¬ 
imate g n \ii=Q obtained using conventional scale setting 
cannot be a RG-invariant. Here, we keep /i arbitrary; its 
value only needs to be large enough to ensure a pQCD 
calculation. The key idea of seBLM is to use the rela¬ 
tion Pi ~ Pq +1 in order to rearrange all the terms at the 


same order following the equivalent /3o-powers, and then 
to eliminate the {/3j}-terms sequentially. 

The first step is to set the scale Q 4 at the NLO level, 
which is determined by eliminating all the {/3i}-terms 
with highest equivalent /3o-P°wer at each perturbative or¬ 
der; i.e. by absorbing all the underlined terms of Ea. (fl9l) 
into the running coupling: 
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g'n = r 0 +ria(Qi ) (l + r 2 [0]a(Qi) + ( /3 0 ■ r 3 [l] + r 3 [0])a 2 (Qi) + ( /?p • r 4 [ 2] + /3i ■ r 4 [0, 1] + ftp ■ r 4 [l] + r 4 [0])a 3 (Qi) + • • 

. _( 20 } 


where the tilde symbol means its value differs from the 
untilded one. The seBLM scale Q \ satisfies 

u 2 

In = ^i,o + Ai,i(/3o ■ a(Qi)) + Ai, 2 (/3o ■ a(Q i)) 2 + ■ ■ ■ , 

Wi 

( 21 ) 

where Aqo, Aiq and A 4i2 are used to eliminate the un¬ 
derlined terms of Eo. ([T9l) at the N 2 LO, N 3 LO and N 4 LO 


levels, respectively. 


The second step is to set the scale Q 2 at the N 2 LO 
level, which is determined by absorbing the doubly- 
underlined terms into the running coupling 


J 


Qn= r o + fia(Qi) ( 1 + a(Q 2 ) ( r 2 [0] + r 3 [0]a(Q 2 ) + (/3 0 • r 4 [l] + r 4 [0])a 2 (<3 2 ) 


( 22 ) 


The seBLM scale Q 2 satisfies 

In yS = A 2 ,o + A 2 p(/3oa(Q 2 )) + 

Qi 


(23) 


where A 2j o and A 2 j are used to eliminate the doubly- 
underlined terms of Eg. GUI) at the N 3 LO and N 4 LO lev¬ 


els, respectively. 


The third step is to set the scale Q 3 at the N 3 LO level, 
which is derived by absorbing the triply-underlined terms 
of Eq. (1^1) into the running coupling 


Qn = r 0 +r 1 a(Qi) (l + a(Q 2 ) (r 2 [0] + r 3 [0]a(Q 3 ) +r 4 [0]a 2 (Q 3 ) H-)) . (24) 


Ea. (l24l) is the final seBLM predictions for the pQCD ap- 
proximant g n , in which the seBLM scale Q 3 satisfies 

In = A 3 ,o H-,, (25) 

b ; 3 

where A 3 ,o is used to eliminate the triply-underlined 
terms of Eo. (l22l) at the N 4 LO level. As with the PMC, 
residual scale dependence remains due to unknown high- 
order {/3i}-terms. 

For definiteness, we will also adopt Q 4 = Q 3 to per¬ 
form the seBLM predictions. The expressions for all the 
A ij coefficients can be found in Appendix [B] In the 
above derivation, we have adopted the equivalent /3q- 
powers [23| to deal with the {/3j}-series at each order. 
We have checked that one can directly replace all /3,; by 
c iP o +1 sca l e setting, and obtain exactly the 

same seBLM predictions. 


3. MseBLM 

As already discussed, the seBLM has some weak points 
which constrain its applicability. In particular, all nj- 
terms are eliminated by the seBLM method without dis¬ 
tinguishing whether those n/-terms pertain to the RG 
/3-function. The seBLM procedure has been illustrated 
by comparing the conformal symmetry limit of QED and 
QCD for the case of the Adler U-function. That is, by 
taking the Abelian N c —> 0 limit of the SU (N c ) group pa¬ 
rameters, such as Cf = 1, Ca = 0, Tp = 1, f abc = 1 and 
d“ bc = l ; together with the condition n/ = 0, the seBLM 
gets the same perturbative series for the D-function as 
that of “quenched QED” (QED with n/ = 0V. Since the 
quenched QED leads to a conformal series j44j , it is thus 
concluded that the seBLM pQCD series is also conformal 
in this case [43[. However, this argument, being based on 
the Abelian limit, is not valid in general; in particular, 
it is invalid for high-order pQCD predictions, because by 
taking the Abelian limit of the QCD series, all the three- 
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gluon and four-gluon couplings are absent [44[ many of 
which also contain conformal contributions. 

Let us discuss these issues based on the RG point- 
of-view: RG-invariance states that a physical quantity 
must be independent of the renormalization scale and 
scheme [451 - 1481 ] . In general, an anomalous dimension 
must also be introduced to ensure the RG-invariance [28l — 
jUj, such as in the case of the Adler D- function. As 
mentioned in the introduction, RG-invariance is bro¬ 
ken at fixed order, leading to well-known residual renor¬ 
malization scale- and scheme-ambiguities. If one re¬ 
quires a fixed-order prediction to satisfy RG-invariance, 
as suggested by Stevenson [ifl (aO ] (called local RG- 
invariance 0 ), one can derive an “optimal scale” and 
even an “optimal scheme” of a process by using the ex¬ 
tended RG-equations. This is the method of the “Prin¬ 
cipal of Minimum Sensitivity (PMS). However, since the 
standard (global) RG-invariance is broken, the PMS pre¬ 
dictions do not satisfy basic RG properties [20], and its 
pQCD convergence is accidental. In addition, it fails to 
achieve the correct prediction of higher-order contribu¬ 
tions when one only knows the NLO correction [54|. This 
limitation also explains why the predicted PMS scale for 
the NLO three-jet production via e + e _ -annihilation does 
not yield the correct physical behavior for low e + e _ col¬ 
lision energy [H2, [53}. A detailed comparison of the PMS 
and the PMC can be found in Refs. [541 ] . 

We find that not distinguishing the n/-terms is in 
some sense equivalent to using local RG-invariance to set 
the renormalization scales. More explicitly, Eqs.(2.12- 
2.16,2.22) in Ref. [HJ agrees with the PMS scale equa¬ 
tion. Thus seBLM will in principle meet the same 
problems of PMS. By eliminating the n/-terms in the 
anomalous dimension function simultaneously with the 
n/-terms for renormalizing the running coupling, the se¬ 
BLM may achieve effective scale of the process for im¬ 
proving pQCD convergence, but it cannot determine the 


correct behavior of the running coupling. 

On the other hand, the /3-pattern and correct {/?*}- 
coefficients using the RG-equation and the standard RG- 
invariance are correctly determined using the PMC. A 
central goal of seBLM is to improve the pQCD conver¬ 
gence; we can retain this goal by a modification which 
will improve its applicability; i.e., we can use the PMC 
method to determine the {/3,}-coefficients to replace the 
seBLM method, while keeping the seBLM procedures 
for eliminating the /3-terms. We will call this modified 
method MseBLM. The MseBLM inherits the main se¬ 
BLM properties, but it avoids the introduction of extra 
MSSM gluinos, thus making it applicable to any process 
and to any order. 

In distinction to seBLM, the MseBLM takes the same 
/3-pattern and the same {/3j}-coefficients as those of 
PMC. It should however only absorb those n/-terms 
which correctly determine the behavior into the coupling 
constant; thus the coefficients for the {/3i}-terms are dif¬ 
ferent from the seBLM ones. We shall show that the 
correct {/3j}-terms are not only needed for achieving the 
optimal running coupling, but they are also important 
for improving pQCD convergence. This partly explains 
why the seBLM and PMC predictions listed in Ref. [23} 
behave quite differently 1 . 


E. Formulas for the Adler D-Function 

It should be emphasized that the anomalous dimension 
function which appears in the Adler D-function has no 
relation to the a s -renormalization of the process; it needs 
to be separately kept fixed during PMC scale-setting in 
order to obtain the correct a s -running behavior 0 - The 
{/3i}-terms that should be absorbed into the running cou¬ 
pling in the Adler D-function can be written as 


D{a) 

D NS (a) 

D s (a) 


1 + 4a + (l2 7 * s + 3/3 0 nf s ) a 2 + (48 7 ^ s + 3/3inf s + 24/3 0 n!f) « 3 
+ (192 7 ? s + 3 / 3 2 nN s + 24/3in^ s + 144/3 0 n^ s ) a 4 
48 7 3 a 3 + (I92 7 f + 144A,nf) a 4 + • • • , 


127r 2 ( 7 (a) - /3(a)—n(a) 


(26) 


(27) 

(28) 


where < 7 / stands for the electric charge of the active fla¬ 
vors. The subscript “NS” and “S” denote the non-singlet 
and singlet parts, da = N c in the fundamental repre¬ 
sentation of the SU(N C ) group. 7 - N ' )S and Hj N)S can be 
found in Ref. 0]. There are no /3 2 -terms at the three-loop 
level, but they are present in the anomalous dimension 


7 ^ s and 7 f and are unrelated to a s -renormalization 2 . 


2 For the anomalous dimension itself, one may need to apply the 

PMC to achieve a better pQCD prediction, which however is out 

of the scope of the present paper. 
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The /3-pattern is different from the seBLM one [23|, in 
which the {^i}-terms from the anomalous dimensions are 
incorrectly included in order to obtain the /3-pattern. In 
the following, we shall show that the correct /3-pattern 
together with correct {/3i}-term coefficients are essential 
for the correct pQCD prediction. 


F. Formulas for R e + e - 


The ratio of e + e annihilation into hadron over muon- 
pairs R e + e - provides one of the most precise tests of 
pQCD. The measured observable R(Q) is defined as: 


Re+e-(Q) 


<r(e + e — > hadrons) 
a(e + e~ —> /i + /x“) 


3 £ 9/(1 + R(Q)), 

f 


(29) 

where Q = yCs = Eqm stands for the e + e -collision 
energy at which it is measured. The timelike i? e+e _ is 
related to the Adler D-function through the equation 


Re+e— (s) 




(30) 


The pQCD approximation for the i?-ratio R(Q) is 


RniQ,^) = ^2Ci(Q,n)a l {iJ,), (31) 

i=l 


where a = a s /4:Tr. Using the Adler D-function, we can 
obtain the four-loop pQCD approximant i? 4 (<3), with the 
initial coefficients of the explicit PMC or MseBLM /3- 
pattern at the scale n = Q given by [56j 


Ci = 3 7l NS = 4, 

C 2 = 12 72 NS + 3/3 0 nN s , 


/ 


c 3 = 48 


7 3 NS 


(E/«) 


2 


3^/9/ 
+3^nf - (7r/3o) 2 7l NS 


V 


and 


/ 


C 4 = 192 


7 4 NS 


+24/3in^ s + 144 


(E s is) 


3 S/ 1 / 

/ 


2 7l 


3/3 2 nf s 


n 


NS 




3 S/9/ 


(32) 

(33) 


24/3 0 nf 


m 


(34) 


-12 (np o y^ - ^ 2 /3 0 /3 l7 r - 3^/3 0 J n^(35) 


The coefficients Ci at other scales can be obtained via the 
scale displacement relation ([5]). 

Following the standard procedures for seBLM, Mse¬ 
BLM and PMC, one can achieve the predictions under 


various scale-setting approaches. The resulting pQCD 
approximation for the i?-ratio is 

n 

Rn(Q,ri = (Q)« l (Qf (Q,/r)), (36) 

2=1 

where Cf s (Q) are free of {/3j}-terms and free of initial 
scale choice. Here the Q f s are the effective scales for 
each scale-setting approach where SS stands for seBLM, 
MseBLM and PMC, respectively. 


III. NUMERICAL RESULTS FOR R e+e - 

In order to obtain numerical predictions for R n , the 
QCD parameters will be fixed using a s (Mz) = 0.1185 ± 
0.0006 [E3l- To be consistent, we shall adopt the nth-loop 
a s -running to do the calculation, and the A^jg is deter¬ 
mined by using the nth-loop a s -running determined from 
the RG-equation. For example, we obtain AL_t = 214 
MeV for i ? 4 by using the four-loop a s -running. It is 
found that the residual scale dependence for all scale set¬ 
tings are highly suppressed; if not specially stated, we 
will take the initial scale fi = Q. For definiteness, we will 
set Q = 31.6 GeV |58| to do our analysis. 


A. A comparison of three-loop R3 



Ci 

c 2 

C 3 

Conv. 

4 

22.548 

-819.496 

PMC 

4 

29.444 

-64.248 

seBLM 

4 

1.333 

-2.298 x 10 3 


TABLE I: Final expansion coefficients C n for R 3 after ap¬ 
plying conventional (Conv.), PMC and seBLM scale settings. 
The T=three-loop a s -running is adopted. Q = 31.6 GeV. 

The seBLM stops at the N 2 LO QCD corrections. After 
applying PMC or seBLM, we can fix two PMC or seBLM 
scales from the known {/3j}-terms, which are 

qPmc = 40 84 GeVj qPmc = 29.79 GeV, (37) 
q seBLM = 24 3 g GeVj QseBLM = 2 52 X 10 " 35 Ge(& 8 ) 

It is seen that the PMC scales are physically reason¬ 
ably, whereas the seBLM scale Q| eBLM is far from the 
pQCD domain, making the final pQCD prediction unre¬ 
liable. This was already observed by Refs. jS, 23], which 
also indicates the importance of the correct determina¬ 
tion of {/3i}-terms. We present the final coefficients C n 
with n = (1,2,3) after applying conventional (Conv.), 
PMC and seBLM scale settings in Tabic [T| One finds 
C n = C n for conventional scale setting. After apply¬ 
ing PMC, the coefficients C n are the scheme-independent 
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conformal terms. After applying the seBLM, the coeffi¬ 
cients C n are the remaining terms by eliminating all rif- 
terms. As shown by TablelJ after applying the PMC, the 
third-order coefficient C 3 is much smaller than the con¬ 
ventional one; in contrast, after applying the seBLM, the 
third-order coefficient C 3 is even larger than the conven¬ 
tional one. We point out that a comparison of the PMC 
and seBLM coefficients and scales after finishing only the 
first step of the procedures, as was done in Ref. [23|, is in¬ 
consistent. 


R 2 


i?3 


Ra 


Kl 


K2 


«3 


Conv. 0.04785 0.04666 0.04635 7.4% -2.5% -0.7% 
PMC 0.04767 0.04667 0.04637 7.0% -2.1% -0.6% 
MseBLM 0.04767 0.04654 0.04640 7.0% -2.4% -0.3% 

TABLE III: Numerical results for R n and K n up to four-loop 
level under conventional scale setting (Conv.), PMC and Mse¬ 
BLM. The value of R\ =0.04454 is the same for all scale set¬ 
tings. Q = 31.6 GeV. 


LO 


NLO N 2 LO total 


Conv. 0.04497 0.00285 -0.00116 0.04666 
PMC 0.04296 0.00380 -0.00009 0.04667 
seBLM 0.04721 

a-BLM 0.04622 0.00018 0 0.04640 

TABLE II: Each QCD loop’s contribution to R 3 under con¬ 
ventional (Conv.), PMC and seBLM scale settings. The total- 
column stands for the sum of all those loop corrections. The 
xBLM results are also presented. Three-loop a s -running is 
adopted. Q = 31.6 GeV. 

We present the separate contributions from the one- 
loop (LO), two-loop (NLO) or three-loop (N 2 LO) QCD 
contribution to R 3 under PMC, seBLM and conventional 
scale settings in Table [H] Due to the elimination of di¬ 
vergent renormalon terms, the PMC pQCD convergence 
is better than that of conventional scale setting, i.e. 

1-^3,pmcI ^ I^pmcI ^ I-^pmcI- 
The seBLM is designed to improve the pQCD conver¬ 
gence by applying the “large /3o-approximation”. How¬ 
ever, the seBLM pQCD convergence for R e + e - becomes 
even worse than the conventional one, and - since it in¬ 
volves an unreasonable small scale Q| eBLM - the seBLM 
pQCD prediction is questionable. To cure this problem, 
a seBLM alteration, i.e. the ai-BLM Q or equivalently 
the one-scale seBLM [23|, has been suggested. Such an 
alteration requires an overall modification of the total 
pQCD prediction by directly requiring the N 2 LO-term 


C 3 to be zero; i.e. Cf 


iBLM 


= 4, C. 


iBLM 


= 1.333 and 


/^ieBLM 
u 3 


= 0. This treatment makes the pQCD prediction 
more reliable as shown by Table [II] However, it breaks the 
expected pQCD convergence since the high-order terms 
C^l M are in general nonzero. Moreover, it cannot im¬ 
prove the seBLM applicability due to the introducing of 
auxiliary fields, which also stops at the N 2 LO level. 

In the following, we shall adopt MseBLM, as an alter¬ 
ation of seBLM, to do a four-loop estimation. 


loop corrections are included, we define a ratio: 

Rn-\-l Rn 
K « — 5 > 


(39) 


where n = (1,2,3). This ratio indicates how the ‘known’ 
lower-order estimate is affected by a ‘newly’ available 
higher-order correction. At the one-loop level, we have 
no information to set the scale for R\, so we take its 
scale as Q and we obtain R\{Q) = 0.04454 for all scale 
settings. Table Hill shows that one can achieve acceptable 
pQCD predictions with increasing loop corrections from 
all scale settings. The K n values for all scale settings be¬ 
have very similarly. The ratio | AC 3 1 for each scale setting is 
smaller than 1 % up to the four-loop level, indicating that 
the four-loop pQCD predictions for R(Q) are sufficiently 
precise. 

LO NLO N 2 LO N 3 LO total 
Conv. 0.04499 0.00285 -0.00117 -0.00033 0.04634 
PMC 0.04290 0.00352 -0.00004 -0.00002 0.04636 
MseBLM 0.04294 0.00352 -0.00004 -0.00001 0.04641 

TABLE IV: The LO, NLO, N 2 LO and N 3 LO loop contribu¬ 
tions for the approximant Ra assuming conventional scale set¬ 
ting (Conv.), PMC, seBLM and MseBLM. The fofabcolumn 
stands for the sum of all of those loop corrections. Four-loop 
a^-running is adopted. Q = 31.6 GeV. 



Ci 

C 2 C 3 C 4 

PMC 

4 

29.444 -64.248 -2.813 x 10 3 

MseBLM 

4 

29.444 -64.248 -2.813 x 10 3 

Conv. 

4 

22.548 -819.496 -2.059 x 10 4 


TABLE V: Final expansion coefficients C n for Ra after ap¬ 
plying the PMC and MseBLM. The expansion coefficients for 
conventional scale setting are also presented as a comparison. 
Four-loop a s -running is adopted. Q = 31.6 GeV. 


B. A comparison of four-loop Ra 

We first present an overview of how QCD loop cor¬ 
rections affect the pQCD estimates. Numerical results 
for R n (n < 4) under conventional scale setting (Conv.), 
PMC and MseBLM are presented in Table Hill To show 
how the theoretical prediction changes as more-and-more 


Next, we present a comparison of pQCD convergence 
assuming various scale settings in Table IIVI The stan¬ 
dard pQCD convergence under conventional scale setting 
is due to the a s -power suppression alone. The PMC and 
MseBLM pQCD series follow the pattern of the stan¬ 
dard pQCD series, but are more convergent. As shown 
in TablefV] after applying PMC and MseBLM, the diver¬ 
gent renormalon terms are eliminated and we get much 
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smaller expansion coefficients C n at higher orders. The 
MseBLM and PMC have the same coefficients C n . This 
explains why a more convergent pQCD series can be 
achieved for both PMC and MseBLM. 



Qi 

Qi 

Qi 

PMC 

41.23 GeV 

36.91 GeV 

171.43 GeV 

MseBLM 

41.03 GeV 

33.61 GeV 

7.23 TeV 


TABLE VI: The determined PMC and MseBLM scales for 
f? 4 . Four-loop a s -running is adopted. Q = 31.6 GeV. 

It is important to find the correct /3-pattern of a pro¬ 
cess. The PMC respects RG-invariance and improves the 
perturbative series by absorbing all /3-terms governed by 
the RG-equation into the running coupling. Thus, the 
PMC eliminates the factorial growth of the renormalon 
nonconformal series and determines and thus yields more 
precise pQCD predictions. The MseBLM adopts the 
same /3-pattern as that of PMC, and as shown by Tables 
m and El its predictions are close to those of PMC. 
There are slight differences for the pQCD series due to 
different way of absorbing the {ft}-terms into the run¬ 
ning coupling; i.e. in contrast to the PMC, the MseBLM 
absorbs the {ft}-terms via the “large /3o-approximation”. 

We have presented the PMC and MseBLM scales for 
i ?4 in Table fVTl one finds different effective running cou¬ 
plings at each order. Following the procedures in Sec. 
II, the PMC and MseBLM scales are themselves given 
as perturbative series. The differences between the PMC 
and MseBLM scales are suppressed by the accuracy of 
the approximation ft ~ /3 q +1 and a s -suppression. The 
scale differences are formally suppressed by the inverse of 
equivalent /3o-powers. In the case of the PMC and Mse¬ 
BLM scales with the highest equivalent /3o-power, we find 
a smaller scale difference. This qualitatively explains why 
the PMC and MseBLM scale differences become larger at 
higher-orders. In the case of the LO scale Q i, it has the 
highest equivalent /3o-power and thus the scale difference 
is the smallest. More explicitly, we have 

A Q x < A Q 2 < AQ 3) (40) 

where A Q t = |Qf MC - Qf seBLM |. 

Finally, we compare the predictive power of PMC 
and MseBLM. We adopt the conservative procedure 
suggested by Ref.jU to predict the “unknown” high- 
order pQCD corrections; i.e., we identify the pertur¬ 
bative uncertainty with the last known order. Thus 
the “unknown” high-order pQCD correction is taken as 

^±|C„a” +1 |MAx) for R n , where \C n a n+1 \ is calculated by 
varying fi £ [Q/ 2,2 Q\, and the symbol “MAX” stands for 
the maximum \C n a n+1 \ within this scale region. The er¬ 
ror estimates for conventional, PMC and MseBLM scale 
settings are displayed in FIG. ([2]). It shows that both the 
PMC and MseBLM errors are smaller than those assum¬ 
ing conventional scale setting. The PMC and MseBLM 
errors tend to shrink much more rapidly with the incre¬ 



FIG. 2: Results for R n (n = 2, 3,4) together with their errors 
(±\C n a n+1 \ MAxj at Q = 31.6 GeV. The big dots, the squares 
and the diamonds are for conventional scale setting (Conv.), 
PMC and MseBLM, respectively. 


ment of pQCD order, consistent with the pQCD conver¬ 
gence shown by Table El 


IV. SUMMARY 

In this paper, we have presented a detailed comparison 
of PMC and seBLM scale settings. Both of these methods 
are BLM-like approaches with the purpose of extending 
BLM up to all orders. 

The PMC provides the underlying principle for BLM, 
which sets the optimal scale for any QCD processes up 
to all orders via a systematic and process-independent 
way. Following the superposition of the RG-displacement 
equation, the PMC also provides a direct explanation for 
the seBLM /3-pattern. The {ft}-terms are then elimi¬ 
nated by the PMC by absorbing them into the renormal¬ 
ization scales consistent with the RG equation. 

Predictions for physical observables should be indepen¬ 
dent of the choice of the schemes or theoretical conven¬ 
tions. The PMC respects RG-invariance, and the final 
expression is scheme and scale independent at each fi¬ 
nite order. There is a small residual scale dependence 
due to unknown higher-order {ft}-terms which is highly 
suppressed even for low-order predictions. The PMC im¬ 
proves pQCD convergence due to elimination of divergent 
renormalon terms. The PMC can also be applied to pro¬ 
cesses with multiple physical scales; e.g., T(IS') leptonic 
decays [jjj. The effective number of flavors n/ is set cor¬ 
rectly at each order of perturbation theory. The PMC 
is consistent with the standard scale-setting procedure of 
Gell Mann and Low in the Abelian limit. 

The seBLM is designed to improve the pQCD conver¬ 
gence, but it cannot be regarded as a solution for solving 
renormalization scheme-and-scale ambiguities. All of the 
seBLM scales are consistent with a (/3 0 • a)-power series, 
and the “large /3o-approximation”. It sets the /3-pattern 
approximately via the relation, ft ~ /3 q +1 , and the {ft}- 
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coefficients are fixed by introducing extra MSSM gluinos. 
The predictions of the seBLM method strongly depend 
on the knowledge of h s -dependent pQCD series and fig- 
dependent /3-functions. At present, the seBLM is only 
applicable to Adler D-function involved processes, and it 
can only be applied to the three-loop level. However, In 
contradiction to its main goal, we find that the pQCD 
convergence of seBLM is questionable, as shown by its 
application to R e + e ~. 

In order to cure the seBLM shortcomings, we have sug¬ 
gested a modification, called MseBLM, by borrowing the 
PMC idea to set the {/3i}-coefficients, while keeping the 
“large /3o approximation” to deal with the /3-series. It 
does not require the introduction of extra MSSM gluinos 
into pQCD calculations. The MseBLM inherits the se¬ 
BLM properties and makes it applicable to any order. 
By taking the four-loop calculation of R e + e - as an ex¬ 
ample, we have shown that the MseBLM and PMC pre¬ 
dictions are numerically consistent; thus more precise 
pQCD approximation and a more convergent pQCD se¬ 
ries for R e + e - can be achieved for both the PMC and the 
MseBLM. This emphasizes the importance of the correct 
knowledge of the /3-pattern and the {/3j}-coefficients. 

In conclusion, the modified seBLM - the MseBLM 
provides a practical approach for improving the 
convergence of pQCD predictions. When more QCD 
loop terms are considered, it can achieve precise predic¬ 
tions consistent with those of the PMC. However, the 
PMC has a rigorous theoretical foundation, satisfying 
all self-consistency conditions from RG-invariance. It 
thus eliminates an unnecessary systematic error for 
high precision pQCD predictions, and it can be widely 
applied to high-energy processes. 
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up to three-loop level. At the scale /z = Q, we have 


ri 


Ml] 
Mo] 
M2] 
Mo, i] 
Mi] 


3 C F = 4 

(Al) 

y - 4 ( 3 ~ 0.691772 

(A2) 

C A C f 1 

(A3) 

3 2 ~ 3 

302 76 

(A4) 

— - y C3 - 3.10345 

^ - 8 C 3 « -1.19979 

(A5) 

^ / 3 80 „ 40 \ 


CA r4 + T C3 ^y C5 ) 



-C F (18 + 52Ca - 80C 5 ) « 55.7005 (A 6 ) 


r 3 [°] = — (523C1 + ^2C a C f - 414C£) 

- . 5/176 \ (£/«/) 

- 72C ^ + s(-- m <dj(^y 


-573.9607- 19.8326 


3 (S/d) 


(A7) 


where Ca = 3 and C F = 4/3 for the SUA3)- group. 
These coefficients can be derived from Ref . 14 il l. There 
is a typo for the coefficient r 3 [l] in Refs. @,01; i.e., the 
term +3/4 there should be corrected to —3/4. The {/3j}- 
coefficients rk[l,m, •••] at other scales can be obtained 
from the above values via the scale displacement relation 
©• 


Appendix A: The {/^-coefficients for the Adler 
D-function 

Following the standard seBLM procedures, we obtain 
the {/3j {-coefficients r^[l, m, ■ ■ ■] for the Adler D-function 


Appendix B: The A ij expressions up to four-loop 
level 
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Ai,o — r 2 [ 1] 

Ai,i = -cir 2 [l] + cir- 3 [0,1] - r 2 [l] 2 + r 3 [2] 

Ai, 2 = i(2c?r 2 [l]-2c?r3[0,l]+3cir 2 [l] 2 -6cir 2 [l]r3[0,l]-2cir 3 [2] 

+2 Cl r 4 [l, 1] - 2c 2 r 2 [l] + 2c 2 r 4 [0, 0,1] + 4r 2 [l] 3 - 6r 2 [l]r 3 [2] + 2r 4 [3]). 

A _ M 1 ] - 2r 2 [0]r 2 [l] 

A 2,0 — -jyvj- 

r 2 [0] 

_ 2cir 2 [0] 2 r 2 [l] - 2cir 2 [0] 2 r 3 [0,1] - Cir 2 [0]r 3 [l] + cir 2 [0]r 4 [0,1] + r 2 [0] 2 r 2 [l] 2 

^2,1 r/~v] 2 

V2 [0\ Z 

-2r- 2 [0] 2 r 3 [2] + r 2 [0]r 2 [l]r 3 [l] + r 2 [0]r 4 [2] - r 3 [l] 2 
r 2 [0] 2 

A _ r 2 [0] r 2 [1] r 3 [0] + r 2 [0] r 4 [1] - 2r 3 [0] r 3 [1] 

AS '° " ’ 


(Bl) 

(B2) 

(B3) 

(B4) 


(B5) 

(B6) 


where Ci is defined via the relation /3j = c»/3 q +1 . 

As an application, we adopt these Acoefficients to 
study the renormalization scale dependence. For con¬ 


venience, we first separate the Q-dependence from the 
/3-coefficient rk[l, m, • • ■ ], where Q stands for the typical 
scale of the process or at which it is measured. 


Mm) = n, (B7) 

Mm) = M0] + /3 0 ^r 2 [l] + In , (B8) 

Mm) = ^3 [0] + /3 0 ^r 3 [l] + 2r 2 [0] In + /3 4 ^r 3 [0,1] + In M j + 0* ^r 3 [2] + 2r 2 [l] In M + In 2 , (B9) 

Mm) = r 4 [0] + /3 q ^3r 2 [0] In 2 M + 3r 3 [l] In M + r 4 [2]^ + /3i ^r 4 [0,1] + 2r 2 [0] In + 

^PoPi ^2r 4 [l, 1] + (4r 2 [l] + 6r 3 [0,1]) In ^ + 5 In 2 + Po (V 4 [l] + 3r 3 [0] In 

( 2 2 2 \ / 2 \ 

r 4 [3] + 3r 3 [2] In M + 3r 2 [l] In 2 M + In 3 + p 2 ( r 4 [0, 0,1] + In MJ . (B10) 


We have implicitly taken rk[l,m, ■ ■ ■] = rfc[Z,m, • • -]| m= q. which are 
For processes with several typical scales, the condition is 
more involved but can be done via a similar way. We can r 2 [l] 

read off the scale-dependent {/3i}-coefficients 77 JZ, m, • • • ], 

r 3 [l] 
r 3 [0,1] 
r 3 [2] 
M 2 ] 

r 4 [0,1] 

Ml,!] 


r 2 [l]+ln^ (Bll) 

r 3 [l] + 2r 2 [0]ln^ (B12) 

r 3 [0,l] + ln^ (B13) 

r 3 [2] + 2r 2 [l]ln^+ln 2 ^ (B14) 

2 2 

3r 2 [0] In 2 M + 3r 3 [l] In M + M 2 IB15) 

r 4 [0,l] + 2r 2 [0]ln^ (B16) 

r 4 [l, 1] + (2r 2 [lj + 3r 3 [0, l])ln ^ 

< B17 > 
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u 2 

r 4 [l] -> r 4 [l] +3r 3 [0]ln^2 (B18) 

„ 2 „2 

r 4 [3] ->■ r 4 [3] + 3r3 [2] In ^ + 3r 2 [1] In 2 ^ 

+ ln 3 ^ (B19) 

u 2 

r 4 [0,0,1] -> r 4 [0,0,l]+ln^ (B20) 

Substituting them into Ajj, we find that, except for 

2 

Ai.o = ^[l] + In q 7 , all other A ij coefficients are free of 


/i-dependence. For Ea. (l2T1) . we obtain 

In ^2 = ^ 2 [1] + Ai ; i(/?o • a(Qi)) + Ai^{Po ■ °(Qi)) 2 - 

*21 

This shows that Q 1 , and thus all the high-order seBLM 
scales, as indicated by Eqs. (I 23I25I) . should be indepen¬ 
dent of the initial choice of scale. This property is en¬ 
sured by the local RG-invariance mentioned in the body 
of the text. 
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